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1. INTRODUCTION 

Use of multipolar electrical installations in the problems of low-depth geophysical observation 
and control allows effective investigating and monitoring the local manifestations of geodynamic processes. 
At the same time, resolution of this method in geodynamic control significantly exceeds the capabilities 
of classical electric prospecting methods with a small number of emitting and receiving electrodes -vertical 
electric sounding and vertical electric profiling [1]. Disadvantage of using a large number of poles 
is the complexity of the measurement process, since each electrode in the course of observations needs to be 
moved separately. Technique of electrotomography (resistivity imaging, electrical resistivity tomography) 
can repeatedly use the same positions ofthe supply and measurement electrodes fixed on the observation 
profile. This leads to a decrease in the total number of operating positions of electrodes with significant 
increase in the density of measurements compared to the conventional method of vertical electrical 
sounding [2, 3]. Currently, there is a wide class of electrotomography devices, manufacturers of which are, 
in particular, Phoenix Geophysics, Iris Instruments, R&D enterprise "ERA", research Institute "GeoTech", 
"North-West". However, for application in the problems of geodynamic control organization their essential 
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adaptation directed to the decrease of influence of natural and technogenic destabilizing factors on 
measurement results is necessary (according to requirements of normative documents of geotechnical 
monitoring organization [4]). 

Developed algorithms for compensation of temperature noise, multiplicative noise of natural 
character and industrial noise in geoelectric control systems can reduce the noise value by 40-60 dB 
and increase the signal-to-noise ratio, which increases the probability of correct detection of negative 
geodynamics or reduces the probability of false alarm. To compensate temperature noise, the efficiency 
of regression processing of daily temperature deviations was proved [5-7]. To reduce the influence 
of meteorological factors of different nature, it is proposed to use automatic calibration of the installation [8] 
and compensation of the geodynamic trend [9]. 

Efficiency of vector measurement systems with registration of phase characteristics of 
electromagnetic field, which are the basis of the phase-metric method of geoelectric control with high noise 
immunity [10], should be highlighted. This is due to the fact that the phase parameters of electromagnetic 
field are less affected by various noise compared to the amplitude parameters of the anomalous component of 
electromagnetic field, usually used to analyze the results of observations [11]. The use of specialized 
software in conjunction with electrotomography equipment allows interpretation of observation data in 
the scope of two-dimensional and three-dimensional models. At the same time, algorithms and programs 
aimed at solving the problems of interpretation of geoelectric control results based on compensation methods 
have not been developed at the moment. In particular, there is no formalized approach to the use of the field 
phase characteristics for interpretation of monitoring data and associated problem of geodynamic 
processes localization. The aim of the work is to substantiate application of the geoelectric control 
compensation method using the field phase characteristics for detection and localization of geodynamic 
processes onthe example of the development of the karst-suffusion cavity. 


2. METHODOLOGY 

Phase-metric recording principle of geoelectrical data of geodynamic control. Work [10] explains 
application of the phase-metric method of geoelectric control, namely using multiple sources as the probing 
signal located in close proximity to the investigated object, and required number of vector measuring sensors 
of the electric field. In the simplest case, two point sources A, B and one measuring sensor O located along 
the AB line and at equal distances from the sources can be used as shown in Figure 1(a). 

When using a quadrature source of the probing field, sources A and B form probing signals shifted 
in phase by 7/2 relative to each other. It should be noted that with a different location of sources relative to 
the sensor, as wellas with multi-pole sensing, phase shifts between the test signals may be different. Each of 
the point sources forms at the point O anelectric field signal of the following form: 


Ey, = Ei, + E Ey = Ei (1) 
Ee =F. + AER,» En - By 


where E” is a normal signal in the absence of heterogeneity; AE’ is an anomalous electric field component 


caused by the presence and geodynamics of heterogeneity. 
As a result of the superposition principle of sources and fields the resulting normalized signal at the 
point O, in accordance with Figure 1(b), has the form: 


H, =(E%./1, + AE,,./I,)sin( at) +E}, /I, + AE,,./I,)sin(at+¢), 


. (2) 
H, =AE,, sin at)/1, +AE,, sin( at + ~)/T, 


A distinctive feature of this method is that the information recorded parameter is not amplitude, 
but phase @, ‘ of the resulting signal, determined in accordance with the following relations at p= 7/2: 


0 
ġ. =arctg Part Abe no , Q, =arctg — 2 (3) 
Ex + AE». AE», 
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Figure 1. Phase-metric principle of geodynamic control 


Geodynamic variations of the studied object are determined by displacement of fictitious sources, 
which leads to an imbalance of the measuring system and registration of the corresponding signal vector in it. 
The phase-metric registration method of geoelectric data is used to monitor geodynamics of near-surface 
inhomogeneities in cases of the need to provide increased sensitivity to special changes in the object of study. 
High efficiency is achieved by increasing sensitivity of the measuring system, initial installation 
and operational positioning of the electrical installation by controlling the sources of probing signals [11, 12]. 

Dynamic control mode of positioning of the phase-metric geoelectric control systems. Algorithms 
for positioning control of the phase-metric measuring geoelectric systems, as well as information processing, 
are based on the assumption that geodynamics of an object can be determined with a sufficient degree of 
accuracy by the results of field registration at N points with coordinates (xi, yi) in the presence of M sensing 
points [11]. In accordance with superposition principle, a vector electric field signal of the following type is 
generated at each i-th measurement point by the source J: 


@.,, =[vi +G,(jo)> BM, expig,), 
k=l 


_ r , 7 , (4) 
yi = V T G, (jo)> BM, exp( JQ; ) 
= 


where Wi, Wi are relative position spatial functions of the probing signal source and measuring sensors; 
Bör» Bir are relative position spatial functions of the imaginary source k, determined by anomalous 


geodynamics of the object under study; /;,@,; are amplitude and phase of the probing signal; Gjj@) is a 
contrast ratio of the imaginary sources. 

The spatial functions in (4) can be determined at the initial system initialization by the method 
of sequential activating of probing signals at the selected poles. To ensure high sensitivity of systemto small 
geodynamic changes of the object, it is necessary to carry out the initial installation and balance it 
on the recorded electric field vectors. This will eliminate from the recorded signals the constant component 
determined by initial inaccurate installation of the probing signal sources and natural misalignment 
of electrolocation signals from the object under study. The balancing condition can be defined 
from the expression ê, = ê, = 0. 

The principle of the system dynamic control is based on the constant adjustment of system due 
to phase control of probing signals. Under condition of initial positioning and that the probing signals 
of multipolar phase-metric systems can be represented with high accuracy either by single harmonic 
oscillations with certain amplitudes, frequencies and phases, or by a set of such oscillations: 


Ae, (j@) =G,(jo)) |) AB; AT, (ja). 


Ae), (jo) = G, (ja)> >A AL (JO) 


j=l k=l 


(5) 


When using a reference quadrature generator in the formation of probing signals of a multipolar installation: 
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AI (jæ) = AL; sin at + AL, cos æt (6) 


where Al eS Al j COS Q, Al i Al š sin 0 are amplitudes of the source j quadrature components. Based 


on (3), we obtain phase geodynamic displacements at the registration points of geoelectric field: 


Re(G,( jo) LAF Aly +Im(G;( jo) DAB A | 


ti ~ RAG Gand Y ARAL, —ImGGand ) ARAL, 
j=l k=l j=l k=l (7) 
Re(G,( a » AL, +1m(G, oae aA i 
Ad, = j=l k=l j=l k=l 


J=l k=l J=l k=l 


Detection and localization of near-surface geodynamic process on the example of karst-suffusion 
cavity. Problem of interpretation of the parameters of geoelectric section based on the observation of several 
parameters of electric field using a multipolar installation can be defined in general as: 


p(r) =argmin ) won +a O° (8) 
A 


where p is resistivity, r = (x y zZ)" is a vector of coordinates of a small fragment of geological medium, A 
is an observed parameter of electric field (potential or field strength, as well as their amplitude or phase of its 
longitudinal or transverse component), wy, is a weight coefficient of the parameter, o, is a mean square error 


(MSE) on the interpretation parameter, a is a regularization parameter; ||€2||? is quadratic norm 
of the regularizing functional containing a priori information about the geoelectric section. MSE of 
interpretation by parameter A is determined by the expression: 


2 
on Iy 


On 5M (A, -A,) (9) 


where M is a number of registration points of geoelectric field, is a weight of the registration poimt of 
geoelectric field with respect to the field source j, and is an exact and approximate value of the field 
parameter at the receiving point, respectively. 

Including weight coefficients is necessary, since in practice it is often difficult or impossible to 
ensure a uniform distribution of electromagnetic field registration points during geodynamic studies. This is 
due to the fact that most geodynamic objects are located in areas of complex development, typical for most 
industrial enterprises or in conditions of limited access. When solving the problems of geodynamic control, 
collection system used is initially stationary and uneven. Due to the fact that distances fromthe object and to 
the source of the probing signal are initially different, registration errors of electromagnetic field at different 
observation points have different variances. Therefore, a distinctive feature of the application of processing 
methods in geodynamic measuring systems, which allows to significantly reduce the error of geodynamic 
estimates, is the use of weight coefficients [13-20]. Weight coefficients are determined based on the 
dependence of the variance of approximation errors of the field parameter on the argument, and accordingly 
determined by localization of the field registration point ri in relation to the probing field source rj 
and the object of study: 


wj =1/ LŒ), SS Wp, = (10) 


n=1 m=1 


Interpretation problems of geological section parameters in general, and localization of 
heterogeneity in particular, are ill-posed, therefore, application of regularization by A. N. Tikhonov is 
necessary [14, 15]. In addition, the problem of lack of a priori information about the section 
for the construction of regularizmg functional Q, characteristic for both electrical prospecting and monitoring 
observations, should be noted. This forces the use of additional simplifying assumptions about the nature 
of the localized heterogeneity. In particular, distribution of the point electrode field in the presence of a karst - 
suffusion cavity is conveniently described by a model of spherical inhomogeneity [21-26]: 
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where is a Legendre polynomial; is a cavity radius; is distance from the sensing point to the center of 
inhomogeneity; is an angle between the direction from the center of inhomogeneity at the point of phase 
control with respect to the direction at the sensing point; coefficient is determined by boundary conditions 
and spatial parameters of subsurface heterogeneity as shown in Figure 2. Also geodynamic changes of cavity 
parameters are designated in Figure 2 as Ad, AO, Ah. 
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Figure 2. Calculation model of karst-suffusion process 


3. SIMULATION RESULTS 


To verify localization possibility of karst-suffusion cavity, simulation of the interpretation procedure 
of the section parameters according to the expression (7) has been done. The model parameters were set: 
installation base L=AB=100 m, probe frequency f=72 Hz, electrical parameters of the host medium 
(moderately moist sand) pl=100 Ohm - m and ¢€1=20, electrical parameters of the karst cavity p2=1000 
o1=100 kOhm - m and ¢2=1, cavity depth h=10 m, cavity radius R=3 m. Distance from the source A to the 
sphere center along the earth's surface was taken to be xl=20 m. Distance x0=AO varied (30 and 50 m), 
which is important for complex construction conditions, when the use of a symmetrical installation 
AO=OB=L/2 is impossible. Amplitudes of the sensing currents were set by expressions [A=] A and to fulfill 
the balance condition of installation in its initial state (in the absence of a karst cavity). 

Figure 3 shows changes in amplitude (Figure 3(a)) and phase (Figure 3(b)) of the anomalous field 
potential component when moving the sphere center along the installation profile from x0 to x=L for 
different values of phase shifts between the probing signals pA B=00; 450: 1350; 900; 1800. Simulation has 
shown a good potential differentiation from the location of heterogeneity. As can be seen from Figure 3(a), 
there is an ambiguity of the potential amplitude on x, which forces the use, along with the amplitude 


component, also the potential phase component. At pAB=00 and 1800 phase differences are manifested only 
in the polarity of the received signal (Figure 3(b)). 


mod, mY 
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Figure 3. Changes in amplitude, (a) Phase, (b) The anomalous field potential component when moving the sphere 
center along the installation profile 
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Figure 4 shows MSE change of the inhomogeneity localization at its displacement along the profile 
with respect to the true position x1=20 m. Localization error is simulated for different weight coefficients of 
the amplitude component (module) of the field w_mod=1; 0.5; O and the phase component 
w_phase=l-w_mod. Best results for both symmetric (Figure 4(a) and 4(b)) and for asymmetric installation 
(Figure 4(c)), characterized by absence of local minima and the smallest global minimum region, obtained 
with simultaneous equal use of the amplitude and phase components w_mod=wphase=0.5. It should be noted 
that for another geoelectric sections or another types of probing installations, the optimal weight coefficients 
may be different. 


oz 
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Figure 4. Localization MSE of the inhomogeneity while displacing along the profile (x1=20 m is the true 
position of inhomogeneity), (a) Symmetrical installation xX0=50 m, pAB=1800, (b) Symmetrical installation 
x0=50 m, pAB=900, (c) Asymmetric installation x0=30 m, pAB=900 


4. CONCLUSION 

Use of vector measurement systems with registration of electromagnetic field phase characteristics 
is the basis of the phase-metric method of geoelectric control, which has increased noise immunity. In this 
paper, application of the compensation method of geoelectric control is done using the phase field 
characteristics for detection and localization of geodynamic processes on the example of the karst -suffusion 
cavity. The principle of dynamic positioning control of phase-metric geoelectric control systems, providing 
increased sensitivity to special changes in the object of study by controlling the sources of probing signals, 
has been described. 

A formalized description of the problem of detection and localization of the near-surface 
geodynamic process based on the interpretation of the geoelectric section characteristics obtained by 
observing several electric field parameters using a multipolar installation is developed. Section characteristics 
are determined by the minimum sum of the weighted mean square interpretation error and the regularizing 
functional containing a priori information about the geoelectric section. 

To verify localization possibility of a spherical karst-suffusion cavity, a simulation of the procedure 
for interpretation of the section parameters based on the developed technique was carried out. Changes in 
amplitude and phase of the anomalous field potential component when moving the sphere center along 
the installation profile for different values of phase shifts between the probing signals is studied. Simulation 
has shown a good potential differentiation from the heterogeneity position. The smallest value 
of heterogeneity localization error for both symmetric and asymmetric installation is obtained with 
simultaneous use of amplitude and phase component in interpretation problem, which indicates the efficiency 
of combining amplitude and phase methods of geoelectric control of geodynamic processes. 
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